Quasiparticle spectrum of the cuprate Bi 2 Sr2CaCu 2 08+<5 
Possible connection to the phase diagram 
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We previously introduced [T. Cren et al., Europhys. Lett. 52, 203 (2000)] an energy-dependant 
gap function, A(E), that fits the unusual shape of the quasiparticle (QP) spectrum for both BiSr- 
CaCuO and YBaCuO. A simple anti-resonance in A(E) accounts for the pronounced QP peaks in 
the density of states, at an energy A p , and the dip feature at a higher energy, Edi p . Here we go a 
step further : our gap function is consistent with the (T, p) phase diagram, where p is the carrier 
density. For large QP energies (E 3> A p ), the total spectral gap is A(E) ~ A p + A v , where A v is 
tied to the condensation energy. From the available data, a simple p-dependance of A p and A v is 
found, in particular A v (p) ~ 2.3 ks T c (jp). These two distinct energy scales of the superconducting 
state are interpreted by comparing with the normal and pseudogap states. The various forms of the 
QP density of states, as well as the spectral function A(k,E), are discussed. 



INTRODUCTION 

A striking feature of the conventional superconduct- 
ing (SC) state is the small number of parameters 
needed for its description. With the knowledge of the 
BCS quasiparticle spectrum, — \Je\ + Ak(T) 2 , re- 
vealing a gap Ak(T) at the Fermi level (efe=0), the 
Hamiltonian is basically known and the magnetic, 
thermodynamic and transport properties can be de- 
rived 0. The gap Ak(T), which vanishes at T c , is re- 
lated to the pairing interaction Vk.k' via the BCS self- 
consistency relation, giving the ratio 2A(0)/fcsT c = 
3.52, in the weak-coupling isotropic case. It is thus a 
scalar order parameter of the transition, a fact that 
has been verified to a high precision |2, y, |j] . 

In the case of high-T c , the probing of the QP spec- 
trum has not led to a solution. Still, a wealth of infor- 
mation on the magnetic field, temperature and doping 
dependence of the SC state has been obtained |5j, |g] • 
The QP spectral function, as probed using photoe- 
mission (ARPES), or the density of states (DOS) as 
obtained by Scanning Tunneling Spectroscopy (STS), 
reveal additional singularities which are at odds with 
a simple BCS d-wave spectrum 0, S H El El El El 
In the DOS (Fig. 1, curve 1), the QP peaks (P) are 
very pronounced and are followed by a dip feature 
(P/) at higher energy (at E = Edip)- Although the 
origin is still debated, there is some strong-coupling 
effect on the q uasip articles : a self-energy is implied 
[3 El ES El El A mean-field approach is insuf- 
ficient in the context of correlated electrons jl9l |2(J , 
coupled spin-charge degrees of freedom [2lt \2'A |23| . 
phase fluctuations [24I Esl Er| , or a competing order 
HIS SHI]- Th us the SC state can no longer de- 
pend on one parameter. 

The main question addressed in this work is can the 
QP spectrum still be described in simple terms (eg. in 
an extended BCS way) and if so, how does it reflect 



the order parameter? To answer, the details of the 
quasiparticle DOS must be understood. As we showed 
in ref . E3 1 the particular shape of the measured spec- 
trum, illustrated by curve 1, cannot be obtained from 
a simple mean-field gap, giving curve 3. However, 
our resonant gap function A(P), described further in 
Sec. I, nicely fits the vari ety of sp ectra published since 

ElEEllllHlSlilllllHElllllialil Despite a 
number of analyses of the QP spectra, taking into ac- 
count the coupling to a collective mode 0, |H 0, |51 
El El S El, the effect of van Hove singularities or 
the particle- hole asymmetry [33ll45ll4ril47j the effects 
of disorder or phase fluctuations \H& |2Tj , the spectrum 
1 is difficult to derive. 
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Fig. 1 : Curve 1 - Quasiparticle DOS as observed by tun- 
neling showing pronounced peaks at E — ±A P followed by 
dips. Curve 2 - Pseudogap type spectrum observed in the 
vortex core, or for T > T c . Curve 3 - Extended BCS (d- 
wave) DOS, but with a larger gap. All three are generated 
by our gap function, Eq. (10). 

The physical parameters of a self-energy, or equiv- 
alent gap function, have yet to be connected to the 
phase diagram. Important leads on the spectral func- 
tion (or DOS), as a function of doping and tempera- 
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ture, have been obtained using ARPES and tunneling. 
First, a well-defined QP peak at the position A p (our 
notation) develops in the SC state and is of d-wave 



symmetry 



develops m tne bl 

in ei m 



A k oc A p cos(k x - k y ), 



yielding the characteristic V shape in curves 1-3. 
However, the ratio 2A p /ItBT c widely departs from the 
BCS value; A p decreases roughly linearly with p, with 
a large negative slope, from underdoped to overdoped 
sides of the phase diagram [6J, |9j, |49j . As is well es- 
tablished, T c (jp) is dome shaped with a maximum at 
Po — -16, where T c ~ 95 K, for P^Si^CaC^Os+s. 
Thus, the position of the QP peak (at E = A p ) is not 
the energy scale the global SC order parameter. 

The QP spectrum has been probed as a function of 
rising temperature |9l|49j,|5JJ,|5jJ,|5j,|5j,|5^|. Above T c , 
instead of revealing the normal state with a metallic 
DOS, the spectrum displays a pronounced but peak- 
less gap of width 2Apc at the Fermi level. This pseu- 
dogap, see curve 2, disappears at the higher temper- 
ature T* , and has possibly the same angular depen- 
dence as Ak- Moreover, A p and A pa have approx- 
imately the same magnitude 0, |HH 0| ■ One finds 
that A p ss 3fcsT*, so T*(p) follows the identical trend 
as Ap in the phase diagram. The challenge is to under- 
stand three contiguous phases (superconducting, pseu- 
dogap and normal). However, the T*(p) curve on the 
overdoped side, where the data is rare, is the subject 
of hot debate [y, |2(J, |2g, y(| , in particular whether or 
not it crosses the T c (p) dome. The T* inferred from 
NMR Knight shift, resistivity and specific- heat mea- 
surements may correspond to still a higher tempera- 
ture, such as the onset of anti-ferromagnetic fluctua- 
tions, distinct from the vanishing of Apq as observed 
by ARPES H3. 

The origin of the pseudogap is a key question about 
high-T c , and is still highly controversial. Theories fall 
into two qualitative categories : Apq is either a pre- 
cursor to, or competes with, the SC state [SJ. The 
first j3 El E3 is immediately compelling due to the 
phase diagram : below T* incoherent pairs are formed 
(thus the pseudogap) and condense at T c . It also fol- 
lows many aspects of the SC state (low carrier density, 
strong coupling energy, small coherence length, sensi- 
tivity to disorder, etc.) and is theoretically tractable 
[24l Ififil I57I Isij l . In the second category, the competing 
order (charge or spin density wave, Varma currents, 
RVB,...[ll|HlIl|29|]) sends T* down, possibly across 
the T c dome, to a quantum critical point above which 
a new phase is formed. Deciding between these two 
categories would be a significant advance. 

Conductance mapping using STS has provided valu- 
able information. The 'normal' state found within the 
vortex core 0, reveals a pseu dogap analogous to 
the one found just above T c [54| . The local vanish- 
ing of the SC order parameter in a vortex core is due 
to the phase singularity; but contrary to the conven- 
tional case, the pseudogap persists. We found a quasi- 
identical pseudogap, in zero magnetic field, caused by 



weak disorder [(Kj. In both cases, A„ and Apq have 
about the same magnitude [13, LMa, IIlA Ml liiil ■ Thus 
an important constraint on the SC gap function is its 
smooth transition to the peak-less pseudogap when 
phase coherence is lost : curve 1 — > curve 2 (Fig. 1). 

Nozieres and Pistolesi have described a super- 
conducting state when a precursor gap, such as Apq, is 
initially present. However, their model implies that the 
SC gap must be larger than the precursor gap. Thus, 
given curve 2 for the pseudogap, curve 3 is expected in 
the superconducting state, not curve 1. Moreover, the 
high spectral weight of the QP peaks (curve 1) leads to 
an apparent paradox: electronic states seem to move 
towards the Fermi level in the transition to the SC 
state 63]. Any complex self-energy or gap function 
must give the correct energy change for the pseudogap 
to SC transition. 

The situation becomes complex in the case of strong 
disorder : local STS mapping has revealed pseudo- 
gap/S C g ap variations at the surface of BSCCO 
|34l I35L l4Ct l61j . Clearly, in this case^there are changes 
of both SC amplitude and phase [2l|, |5jJ |6j|, |6ja, |6j| . 
In Pb-substituted BSCCO the disorder causes 'super- 
conducting' islands to form where the intensity 
of the spectral fine structure correlates with the de- 
gree of long-range order. Recent theoretical work by 
Atkinson [21| and experimental STS by Fang et al. |40( 
corroborate this conclusion. Even if the inhomogene- 
ity is not an intrinsic property 0, H3j E3 , the cuprate 
SC state is sensitive to local perturbations |64ll65ll66| . 
and the attenuation of the spectral fine structure is a 
clear sign. Then, several parameters are needed for the 
interpretation of the tunneling DOS. 

In this Article, we analyze the QP spectrum, E^, as 
inferred from the sharpest tunneling DOS of BSCCO. 
The detailed DOS shape (pronounced QP peaks, fol- 
lowed by dips at higher energy) is due to a single reso- 
nance in our energy-dependant gap function, Section I. 
Such a resonance in the QP spectrum is possibly the 
coupling to a collective mode [I5LH7II31I 

111 ei Elm 

of the same origin as the resonance seen using inelas- 
tic neutron scattering |6S| and which scales with T c : 
fl ~ 5.3 fcs T c . Zasadzinsky et al. studied the dip 
position as a function of doping using strong-coupling 
theory [J^, and suggested that the dip energy is also 
related to T c . Since the origin of the resonant gap func- 
tion (or QP self-energy) is still unknown, our aim is to 
show how the basic parameters depend on the carrier 
density p. In Sec. II we find that one of the energy 
terms, A v , is compatible with an order parameter: it 
is proportional to ks T c {p). The predicted shape of the 
QP density of states is then studied as a function of 
p. We treat the transition to the PG state (Sec. II), 
and the role of the two distinct energy scales (Sec. III). 
Finally, the QP spectral function and self-energy are 
discussed (Sec. IV). 
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I. SUPERCONDUCTING GAP FUNCTION 

Here the quasiparticle DOS, having the essential 
characteristics of the observed STS conductance spec- 
tra, is derived. In Fig. 2, we show such a spectrum 
obtained on BSCCO, near optimally doped, from our 
group (2a) which is c omp ared to (2b), a spectrum from 
Pan, Hudson et al. [12|. One can again see the pro- 
nounced QP peaks, the steep slope on the outer side 
of each peak, followed by the dip feature previously 
described. In this Section, we focus on these main 
aspects of the DOS; the questions of the background 
slope (as in Fig. 2b), the Fermi surface anisotropy, 
the particle-hole asymmetry and the van Hove sin- 
gul arit y , have already been given extensive treatment 
Hi EE IH H IH. The detailed fits of Fig. 2, solid 
lines, give the key parameters of our SC gap function, 
without such considerations. 

Expression for the QP-DOS 

Consider th e sp ectral function, as measured by 
ARPES 0, H El EJ HI HE E3, in the two- 
dimensional model with k — (fc, 6) the wave vector 
in the ab plane : 



A(k, E) = - Im G(k, E) 

7T 



(1) 



where G(k, E) is the single-particle Green's function. 
The superconducting DOS is then [HEIl : 

N S (E)=Y / A(k,E) 

k 

which is measured in the tunneling experiment. Con- 
verting sums to integrals in the usual way : 

N n (0) >' 2:v 



N S (E) 



n M old [°° de k A(k,E) (2) 

'-I* JO J-oo 



where and N n (0) are the normal excitation spec- 
trum and Fermi-level DOS, respectively. This expres- 
sion ignores the effect of the Fermi surface anisotropy 

|46l l69l| . In the case of an ideal quasiparticle with zero 

71 



lifetime broadening 
the propagator is |3. 

G(k,E) = 



and dispersion E k = ye% + 
E3: 



E-Eu 



iO- 



E + E u + iO- 



with Uf., Vk the usual coherence factors. Then, A(k, E) 
is just: A(k, E) = u\ S(E - E k ) + v\ 6(E + E k ) where 
each term corresponds to a quasiparticle added (E > 
0) or removed (E < 0), respectively. 

For example, with fixed E > 0, the integral in 
(2) picks up two poles at ±efc of am plit ude uf,(— £&) 
and ul(+ek). As is well known |sl l7lj. the coher- 
ence factors disappear in the symmetric case, since 
ul(-ek) + u|(+e&) = 1. We shall thenceforth ignore 
them, and consider for both tunneling and ARPES 



that the values ie^ are equivalent. Using Eq. (2) 
the expected result is obtained : 

f27T 



N S (E) =N n (0) 



1 

2^ 



d6 



dE v 



(3) 



E k =E 



It is then convenient to define the partial (one dimen- 
sional) DOS at the angle 6 : 

N S (E) = / n s {E,6) dO (4) 
Jo 

and, in the extended BCS <i-wave model [USE!]' the 
partial DOS is : 

N n (0) E 



n s {E 7 6) = 



2n 



E 2 - A 2 p cos 2 (26») 



(5) 



Expressions (4, 5) can be used to generate the curve 3 
in Fig. 1 which, as we have stressed, fails to match the 
strong spectral weight of the QP peaks, and the dip 
features, obvious in the experiments (Fig. 2). 
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Fig. 2 : Dotted lines - local STS conductance spectra. 
Solid lines - fits using the model (9, 10). 

(a) Symmetrized spectrum from ref . [Tl| . with back- 



ground removed. Here A p = 31.5 and Edip 
the two free parameters are : A ¥ 



64 are fixed; 
18.5 and E Q = 50, all in 
meV. Also, Eo = Ao , n = A v and A=2, as seen in Fig. 3. 
The broadening values are T = .08 A p and S — 0.1. 

(b) Spectrum from Pan et al. |T^|. The gap function (10) 
is used but the parameters, A p (p) and A v (p), are from the 
p-dependent fit to the data (Fig. 4a). The best fit is for 
p = .175, with the background slope added. All parameters 
are nearly identical with (a) except _Eo — 1.06 Ao and T = 
.03 A p . 

In the general case, beyond the mean-field approach, 
one must evaluate : 

N n (0) 



n s (E,9) 



^ fde k A(k,E) 
In J 



(6) 
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assuming a suitable model for A(k, E) El E 
E3- l32l |42| . Our approach is to consider that strong- 
coupling modifies the electron-electron interaction, but 
without the retarding effects that would occur in the 
case of phonon-mediated pairing Thus we write 
Ak — > Ak(-Ek), and the new dispersion law is : 



between the QP peaks and the dip position; an esti- 
mate is Eq ~ (A p + Edi P )/2. This condition ensures 
that the QP peaks are reinforced, where the cxtrcmum 
of —dA/dE is positive, and gives the dip at a higher 
energy, where —dA/dE is a negative extremum (see 
Fig. 3). 



4 + A k (£ k ) 2 + l r 



(7) A(E) 



where we add T, the lifetime broadening introduced by 
Dynes [7j| . Assuming particle-hole symmetry, we use : 



A ^ = k Im E-lk 



iO- 



CS) 



and, performing the complex integration in (6), we ob- 
tain for the partial DOS : 



n s (E,( 



N n (0) Re E- 



iT 



2vr 



^{E-zTf-A^Ef 



0) 



Here Ak, evaluated at the pole (E^ = E), is a func- 
tion of 8. For a constant gap Ak = A p , we get back 
Dynes's formula for the BCS DOS with the lifetime 
broadening. In the general case, (9) is the basic equa- 
tion for the quasiparticle DOS in our approach, once 
integrated over 6. It contains a new term in the nu- 
merator, —A^dA^/dE which is responsible, as already 
shown in f° r the distinct modification of the DOS 
seen in Fi g. 2 . It can be used to match the tunne ling 

data ofJiMIIlIillllHIMIlllllHIMIiilifl 
MMm 

Superconducting gap function 

We now consider the particular gap function ap- 
propriate for fitting the data (solid lines in Fig. 2). 
For the superconducting state, we have made a slight 
change in notation with respect to ref. llj. Assum- 
ing Ak(-E) = cos(29)A(E), the gap function along the 
anti-nodal direction is now written : 



A(E) = A p + A v (l - g(E)) 



(10) 



where A p and A v 



are constant parameters and g(E) 
is a simple Lorentzian : 



9(E) = A 



(E - E ) 



(11) 



having the standard parameters. Thus the second term 
in (10) is an anti-resonance (a local decrease in the pair 
potential). It provokes an additional peak, and dip, in 
the DOS near the two possible extrema of dg(E)/dE 
(see upper panel, Fig. 3). 

Consider the problem of fitting the DOS in some 
systematic way : there are a priori 5 parameters (the 
role of T will be discussed subsequently). We com- 
mence with the data of Fig. 2a, where the spectrum 
is symmetrized and the background removed. From 
our previous work [llj . the resonance energy must lie 
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Fig. 3 : Characteristics of the SC gap function and their 
link to the DOS along the anti-nodal direction (0 — 0). 

Upper panel - Real part of A(E) from the fit of spectrum 
Fig. 2a. Note the anti-resonant shape, with the minimum 
at the energy Eo = Ao, while both amplitude and width are 
2A V . The intersection point, along A(E) = E, gives the 
QP peak at A v . For large energies, the pairing interaction, 
is A(E) ~ Ac, = A p + A v . 

Lower panel - Partial DOS showing a sharp QP peak 
(P), followed by the dip (D) at the higher energy Edip- 
The derivative of the gap function, as in Eq. (9), reinforces 
the quasiparticle peak (negative extremum) and causes the 
dip feature (positive extremum). 

Since A p is defined as the QP peak position, easily 
estimated from the data, we have the further condition 



that: A(A p 



also illustrated in Fig. 3. 
1, or : 



(10) and (11), this gives g(A p 

(E - A p ) 2 



A = l 



Using 



(12) 



as a constraint on the parameters. A second constraint 
is found by writing the dip position, E^ p , using the 
analytical expression of the DOS. Considering the ex- 
tremum of Ag?A I dE leads to the good approximation : 
Edi V — Eo + .73rj. Therefore, taking Edi P and A p to be 
known, the fit to the spectrum of Fig. 2a can be done 
by the variation of only two free parameters : Eq and 
A™. Their final values determine the precise concavity 
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in the DOS, between the QP peak energy and Edip, 
which is a priori unknown. 

For the fit of spectrum 2a, using A p = 3f.5 and 
Edip = 64, we obtain the values Eq = 50, and A v = 
1 8.5, all in meV. The problem is thus simplified by the 
following outcome : Eq = A p + A v = Ao, with 77 = A^ 
and A=2. 

The corresponding function A(E) is depicted in the 
upper panel of Fig. 3, with the partial DOS along 9 = 
for direct comparison (lower panel). One can see that 
A(E) has a minimum at Eq = Ao, where its value 
is A p — A^, then it increases towards the asymptotic 
value A p + A v , when E ^S> Edip- Note that the total 
amplitude of the resonance, 2A V , is identical to its 
width. In short, we obtain : 



A 



Ap = Eq 



Ap = rj = A ¥ 



and we propose that these relations should scale when 
the carrier density, p, varies (Section II). The dip po- 
sition, near to A p + 2A ¥ ,, is more precisely : 



E d 



ip 



Ao 



.73 A u 



Ap + 1.73 A v 



so that E, 



dip 



A p oc A v . The latter parameter thus 
plays a fundamental role in our model. 

Broadening parameters 

Two broadening parameters in the final fits of Fig. 2 
are used. First, the iT introduced by Dynes to treat 
a finite quasiparticle lifetime, artificially displaces the 
pole in the spectral function (8) off the real axis. It 
is equivalent to a Lorentzian broadening of the DOS, 
of full-width ~ 2r, affecting the states at all ener- 
gies. Consequently, in the lower panel of Fig. 3, vir- 
tual states lie within the gap of the partial DOS (for 
E < Ap). As we discuss in a recent paper [tJJ . sev- 
eral factors can contribute to T both intrinsic (inelastic 
scattering, many-body effects,...) and extrinsic (high- 
frequency noise,...). In Fig. 2, the T value was adjusted 
to fit the zero-bias conductance. 

The QP peaks, as calculated using (10), are initially 
higher than those of Fig. 2. Since T is fixed, a second 
broadening parameter is introduced : we replace A p 
by a complex number with a small imaginary part, 
Ap — * Ap(l — iS). This has a major effect on the QP 
peaks, but a small one on the remaining spectrum. 
Intuitively, the imaginary part represents a smearing 
of the value of A p , such as in the case of gap anisotropy 
[Til l or Doppler shifts due to supercurrents [l^ . While 
this ansatz is used to model the pseudogap in Section 
II, it has a negligible influence on the fit parameters 
deduced here. 



II. CONNECTION TO THE PHASE DIAGRAM 

The fit to the BSCCO spectrum, Fig. 2a, leads to 
the simple result that the width, amplitude and po- 
sition of the anti-resonance are all simply related to 



the quantity A^; assuming A p to be known. The QP 
spectrum therefore depends on only two energy scales. 
This section is devoted to their possible link to the 
phase diagram : we infer how the parameters of the 
SC gap function change with the carrier density. 

Doping dependance of the energies 

The foregoing suggests that A v must have a new 
meaning. Consider the asymptotic value of A(E) for 
large E : A(E) ~ A p + A v = A for E > E dlp (Fig. 3). 
A(E) is then constant up to some higher cut-off en- 
ergy, as in BCS theory. Along the anti-nodal direction : 
E^ « yjej, + Aq, and A is thus interpreted as the total 
spectral gap in the SC state, even though the QP peak 
remains at the smaller energy A p . The pseudogap in 
the vortex core 0, EH, EH , where phase coherence is 
lost, takes on a value Apq ~ A p and, aside from the 
thermal broadening, the same holds for the pseudogap 
just above T c [H |Hl H3 . 
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Fig. 4. Doping dependence of the two basic parameters 
of the gap function, A p (p) and A v (p), their sum, Ao(p), 
and the dip position Edi p - 

Upper panel (a) - A p (p) line : Linear fit to the QP peak 
positions taken from ARPES data J|j4^j5|L SIS/SIN tun- 
neling H m 0U E3. and STM ETOlll Ei El E3- 

Edip line : Best quadratic fit to the dip position, as es- 
timated from STM data ETJ, El E El El @ and 
SIS tunneling [jj IT3I E3- The solution of the fit gives 
A v (p) = 2.3 k B T c {p) and E = 1.06 A (p). 

Lower panel (b) - Analogous diagram to (a) but with 
parameters normalized to the gap energy A p : simple linear 
laws are obtained, Eqs. (14-16). The principal energy scale 
of the SC state, and determining the DOS fine-structure, 
is A lp (p). No evidence for a critical point is seen for 0.05 < 
p <0.25. 
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We now suggest that A v is tied to the condensation 
energy. The energy (per pair) of the SC state is then 
~ — Am, with respect to the PG state, but it is ~ — Ao 
with respect to the normal state. As will be discussed 
in Sec. Ill, the integration over energy states, involv- 
ing the full Ak(£'k), is necessary to obtain the precise 
energy changes. 

We can thus put : A v = C ksT c where C is to be 
determined. From the fit of Fig. 2a, using T c ~ 90 K, 
and A v = 18.5 meV, gives C ~ 2.4. 

With this result, one could use scaling arguments 
to infer the behavior of the parameters as a function 
of p. However, in view of the dispersion of the tun- 
neling and ARPES data, using a single spectrum and 
one value of T c is restrictive. We take a different ap- 
proach by first plotting in Fig. 4a a larger set of data of 
the QP peak position, A p . Assuming that the critical 
temperature T c (p) is known [6j, one can deduce the 
linear behavior of A p (p) by a quadratic fit 0. The 
data are from the SIS (SC-vacuum-SC) break junction 
0- EI EH HI] 7 which probes the electronic structure 
deep in the sample; SIN tunneling data, as well as 
ARPES data, are included HEllHi. The dip posi- 
tion (Edip ), a s seen in many different tunneling exper- 
iments EElElElElllilll,E3,E3, is also plotted 
in Fig. 4a. The points evidently follow a continuous 
curve throughout the phase diagram. We propose to 
find this curve via a best quadratic fit using the gap 
function, which fixes all the parameters. 

The two energy terms are thus A p (p) and A v (p) = 
C fcsT c (p) where the best value of C is to be deter- 
mined. In order to calculate Edi P we could use the 
previous condition Eg = Ao, which led to E d i P — 
A p + 1.73 A^. However, the spectrum Fig. 2b, from 
Pan et al., cannot be precisely fit with the resonance 
energy exactly at Ao. This is the essential difference 
between the two spectra (2a and 2b); the latter dip 
position is slightly at a higher energy. 

From our study of YBCO in [ll|, we noticed the 
resonance energy can be larger than Ao . We thus write 
Eq = AAo, in the general case, and determine the value 
of A at the same time as C . The expression for the dip 
energy is then : 



E dip = AA + -73 (AA - A p ) 



(13) 



and all other parameters retain their previous meaning. 
Putting A p (p) and A v {p) = Ck B T c (p) in (13), and fit- 
ting the data of Fig. 4a, yields : C = 2.3, A = 1.06, and 
the continuous curve Edi P (p). The two values differ by 
about 5% from the single-point estimate and, consid- 
ering the uncertainty in the data, the earlier observa- 
tion in 0| that Eq ~ Ao is maintained. More signifi- 
cantly, we find that A v (p) = 2.3fcsT c (p) is consistent 
throughout the phase diagram, from underdoped to 
over doped sides. 

In Fig. 4a the total spectral gap Ao(p), and the am- 
plitude A ip (p), are plotted as a function of p. One 
observes that A^p) merges smoothly with A p (p) on 



the overdoped side. Consequently, A (p) is a smooth 
convex function of p ranging from A p to 2A p . Extrap- 
olating Edi p (p), it varies from ~ A p to ~ 3A p in the 
same range of p. At optimal doping, p = po, the values 
are A v ~ A p /2 and A ~ 3A p /2. The dip position is 
at E d ip — 2A p , in agreement with ref. |l3| . 

Simple trends are found by plotting the parameters 
as ratios with respect to A p , Fig. 4b. We see that 
A v /A p , Ao/Ap and E d i P /A p are increasing linearly 
as a function of p. Moreover, if we write p as the ex- 
cess carrier density from the minimum value of .05, i.e. 
the doping at the SC onset, we obtain : 



Ay 

A p 
Ao 
A„ 



P 

2p 



E, 



dip 



i + 
1.1 + 



v 

2p 
0.9p 



Po 



(14) 
(15) 
(16) 



where po — -11 is the optimal doping. The first two are 
obtained from the Taylor expansion of A v / A p while 
the third is from Eq. (13), using the value A = 1.06. 
Thus Edi p /A p shows nearly twice the slope as A v / A p , 
and varies from about ~ 1 to ~ 3 in the complete 
doping range. 

Such a straightforward relationship between the pa- 
rameters was not expected. If A v is indeed the 'con- 
densation amplitude', it continuously increases relative 
to the precursor gap, A p , throughout the phase dia- 
gram. Since the latter decreases linearly with doping, 
the order parameter is parabolic shaped : 



A v = 



pAp(p) 
2po 



2.3k B T c (p) 



(17) 



With the previous hypotheses, Eq. (17) expresses a new 
precise relation between the QP peak positions and the 
order parameter of the SC transition. 

To conclude the discussion on the results from the 
fits, it is remarkable that the gap function, A(E), as 
displayed in Fig. 3, scales perfectly as a function of p 
through the variation of its amplitude, A M] and the 
resonance energy, AAo. As with refs. [9j,ll3|, we find 
no abrupt change in the QP spectrum, nor in its un- 
derlying parameters, while spanning the carrier con- 
centration. We conclude that a critical point, if there 
is one, is situated at the right end of the T c dome. 

Shape of the quasiparticle DOS 

The gap function A(E) is now uniquely determined 
for the range of carrier concentration of interest (.1 < 
p < .24). It is then possible to fit the QP DOS with es- 
sentially one free parameter (p) , aside from the broad- 
ening and the background slope. The spectrum of 
Fig. 2b, from Pan et al., was fitted by adjusting the 
value of p, the final value being p = 0.175 (see Fig. 5). 
This could be an approximate measure of the local 
value of the doping at the surface of the sample, but 
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the background slope adds some uncertainty. Our ob- 
jective here was mainly to reduce, as far as possible, 
the number of free parameters. 




-60 -30 30 60 90 
Energy (meV) 



Fig. 5. Plot of the quasiparticle DOS evolution as a func- 
tion of p, the carrier density. The values of the parameters 
used in the gap function are fixed from Fig. 4a : the QP 
peaks (A p ) follow a linear trend, while the dip position, 
Edip, lies along the curve given by (16). The difference be- 
tween the dip and the peak positions is ~ 2A ¥> (p). Note 
the detailed DOS shape, varying significantly from under- 
doped to overdoped sides. 

The variation of the QP DOS shape as a function of 
p is shown as a surface plot in Fig. 5. As expected, the 
peak position follows a linear law, while the dip follows 
the curve given by 

£ di p~l.lA p (p) + 1.8A v (p) 

as a direct consequence of the gap function (10). The 
detailed shape of the QP DOS is not the same on each 
side of optimal doping due to the relative position of 
Ao compared to A p . Indeed, on the extreme under- 
doped side, where A p > A v (p), the outer slope of the 
QP peaks is more rounded, with a slight negative con- 
cavity between peak and dip. On the overdoped side, 
A p ~ A v (p), the DOS has sharper peaks and dips, but 
also a positive concavity between them. A closer look 
at the extreme overdoped case reveals a new singular- 
ity, near Ao (a small kink). The kink is due to a new 
extremum of dA(E) 2 /dE, but this prediction remains 
to be verified experimentally. 

The question as to whether we obtain 'BCS be- 
havior' on the overdoped side can be addressed. In 
one sense we do: this limit gives A p (p) ~ A v (p) ~ 
2.3fcsT c and a single parameter is then tied to T c . 
However, the total spectral gap A becomes twice too 
large (~ 2A p ) so that the dip still persists, near ~ 3A p , 
which contradicts the BCS <i-wave case. Thus, the 
precursor gap maintains its influence when p — > 2p . 
In fact, the results of Figs. 4 and 5 clearly show that 
the gap function is non BCS-like throughout the whole 
phase diagram. 

Our model allows an order of magnitude for the 
value of T* , defined by the vanishing of the precur- 
sor gap. Using the data of d,0)H3l> we nn d roughly 



that 2.8 k B T* ~ A p . Combined with Eq. 14, we get : 
T*~1.6(po/p)T c 

assumed valid for T* >T C . It follows that the two tem- 
peratures merge at p = 1.6po — -18, perhaps remain- 
ing merged for p > l.Gpo, but this is an open question. 
Note that the excess doping level at the (T*,T C ) cross- 
ing point is .18, while the absolute level is .23, i.e. well 
into the overdoped region of the phase diagram. 

Model of the pseudogap state 

We propose a phenomenological description of how 
the gap function should evolve when SC coherence is 
lost, i.e. the pseudogap state. Since the model has 
fixed parameters for the SC state, we are allowed one 
additional parameter to 'force' the transition to the 
PG state. There is very little data on the QP spectral 
change, and we try only to be qualitative; no precise 
fit to a PG-type spectrum is done. 








Temperature 

Fig. 6. Extension of the model to finite temperature 
showing the T-dependence of a and (3. 

Recall that the T-dcpendant BCS gap can be writ- 
ten A(T) = /3(T)A(0), where /3(T) is a solution to 
the finite-temperature gap equation and which verifies 
(3(0) — 1 and (3(T C )=0. A numerically exact, or simple 
interpolation formula, can be found for /3(T). In an 
analogous way, we write : 

A(E) = a A PG (l-/3)+f3 A SC (E) (18) 

where Asc(E) is the previous gap function and /3, a 
vary smoothly (but arbitrarily) from 1 to as a func- 
tion of the strength of the perturbation (temperature, 
magnetic field, etc.). We assume a and (3 vanish at 
T* and T c , respectively, see Fig. 6. The three region 
phase diagram is a consequence of these two parame- 
ters, while Ape and Asc are the characteristic zero- 
temperature energies. In such a way, the three states 
are evidently : 

(3 ^ a^0 SC state 
(3 = a^0 PG state 
(3 = a = Normal state 

In the case where A pq is strictly identical to A p of the 
previous section, we get the form : 

A(E) = a A p + (3 A v (l- g(E)) (19) 
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so that the 'condensation amplitude' is (3 A v and the 
'total' spectral gap is A = a A p + (3 A v . 

The QP-DOS shape corresponding to the PG state 
is an unsolved problem, but there are many models 

EM IS3, HE Q IiE S p3, EM 123, EE E3, m, E3, 

V76L |23 ■ There is a stron g in dication of a broad-band 
self-energy in this case [7U431 IH1 IH1 lol. Id^ . Further- 
more, as discussed in the introduction, the PG-DOS 
is characterized by the absence of quasiparticle peaks. 
Atkinson et al. [ill . 16^ have extended the earlier work 
of Huscroft & Scalettar E3 that a non-condensed Bose 
state produces a peak-less pseudogap in the DOS. This 
has recently been applied to the case of spatial disor- 
der |23. Franz & Millis ^E|j m an altogether different 
approach, have found a pseudogap due to spontaneous 
currents, also typical of a fluctuating pair system. 

Our exceedingly simple model is to consider Ape 
with a finite imaginary part : 



A 



PG 



A p (1 - iS') 



(20) 



and to use Eq. (18) for the gap function. With (3 low- 
ered from 1 to 0, the SC quasiparticle DOS evolves 
smoothly into the PG shape (see Fig. 7a). This type 
of DOS modification is seen in the STS experiment as 
the tip is scanned within the vortex core |39l l59j ] . 
It is also seen as the effect of weak disorder |60j . Here 
5' = .2, which has a small, but non-negligible, influence 
on the energy changes previously discussed. Notice in 
Fig. 7a that the dip features disappear as (3 — > 0, along 
with the QP peaks. 
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Fig. 7. Variation of the DOS with j3 varying from 1 (SC 
state) to (PG state). The parameters are p=.18, giving 
A p = 30 meV and A v = 17.8 meV, in Eq. 18, and 5' = .2, 
in Eq. 20. The broadening parameters for the SC state 
are: 8 = .06, Y = .03 A p . In (a) ReA PG = A„, in (b) 
ReApG = 40meV, i.e. larger than A p . 

These variations simulate the loss of SC coherence in the 
case (a) in the vortex core [T21 159^1 and (b) due to strong 
disorder |H EH Ell- 



in Fig. 7b we let (3 vanish but here Re Ape > A p : 
the pseudogap is slightly larger than the quasiparti- 
cle peak energy (40 meV as compared to 30 meV). 
Now the QP peaks move progressively out to a higher 
energy, as expected from Eq. (18), while being atten- 
uated. This particular evolution is strikingly similar 
to STS observations 0, |3f| |3(| , and is generally 
attributed to disorder (inhomogeneity) . In our model, 
the decrease of doping level on the local scale could 
give this effect, the result being a wider pseudogap. It 
corresponds to the progressive lowering of the p-value 
in the phase diagram of Fig. 4a. However, a disorder 
potential is expected to directly affect the SC ampli- 
tude, hence ReApc could vary in different regions of 
the surface EH- Consequently, the SC state to PG 
state transition energy depends on the local change of 
the precursor gap, (Re Ape — A p ), and on A v . An es- 
timate of this energy is given in the following section. 




-100 100 

Energy (meV) 

Fig. 8. Simulation of the tunneling conductance with 
temperature. The value of f3, varying from 1 (SC state) to 
(PG state), was inferred from the data |54{|. shown in the 
upper inset. The parameters are p—.14, giving A p — 42 
and A v = 17.6 meV, in Eq. 18; all others are the same as 
in Fig. 7a. Here we assume ReApc = A p , but a in Eq. 18 
decreases slightly from unity (see inset). The highlighted 
spectrum corresponds to T c . 

The well-known experiment of Renner et al. [3] 
gives the quasiparticle conductance spectrum as a 
function of rising temperature. In Fig. 8 we plot a sim- 
ulation of the temperature dependance of the QP spec- 
tra using Eqs. (4,9) and the gap function (19). The 
standard broadening of Nc(E) due to Fermi statis- 
tics is taken into account 74]. The value of (3{T) is 
inferred from the data (upper inset, Fig. 8) and it van- 
ishes at T c ~ 89 K. The p value is chosen such that 
A p (p) ~ 42 meV. As is well established from many 
experiments, the PG persists well above T c , up to the 
temperature T* previously discussed. Our model, as 
shown in Fig. 8, gives a simple qualitative description 
of the QP DOS with rising temperature. 
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III. INTERPRETATION OF THE ENERGIES 



We turn to the question of interpreting the two en- 
ergy scales (A p , A v ) that are the main ingredients 
of the gap function A(E), as expressed in (19). The 
foregoing work shows that, for large QP energies and 
along the anti-nodal direction, the total spectral gap 
is A(E) ~ A = A p + A Vl where A p defines the QP 
peak positions, and the new parameter, A v , is propor- 
tional to T c (p). We now consider the energy changes 
involved in both the normal to SC and the PG to SC 
transitions. In the latter case we use (20) for the PG 
gap, with the assumption : Re Ape — A p , with a small 
imaginary part. 

Consider first the BCS case of a constant isotropic 
gap of value A in the SC spectrum. The approximate 
change in energy is given by the familiar integral : 



(E N - E SC ) BCS = JV„(0) / — / de k -f - |e 



2n 



cie 



2tt J \E k 
ny ' J 2nJ k 2E k 



(21) 



The first and second terms [78j are the change in ki- 
netic and potential energies, respectively, while the 
variation of the chemical potential is neglected. The 
high energy cut-off is unnecessary since the integrand 
vanishes faster than l/e^. The standard result is: 



(E> 



Esc) BCS = g ^™(0)A 2 



or equivalcntly, putting N p = N n (Q)A/2 as the number 
of pairs, 



(En — Esc) 



BCS 



N, 



(22) 



The conventional interpretation is therefore straight- 
forward : the gap in the QP spectrum is equivalent to 
the energy gain per pair in the N — > SC transition. 
The latter intensive quantity is related to the mean- 
field critical temperature, A = 1.76 k B T^ F . 

In the unconventional case, we must go beyond 
the mean-field due to three possible effects (a) weak 
phase-stiffness and 2D electronic structure (Kosterlitz- 
Thoulcss limit), (b) proximity to Bosc-Einstein con- 
densation and (c) coupling to a collective mode (ad- 
ditional degrees of freedom) . Denoting E v the energy 
gain due to one of these processes : 

(E N - E sc ) = (E N - E SC )' BCS + E v 

where (En — Esc)bcs ^ s gi ven by (21), but using the 
complete gap function A{E k ) (i.e. an extended BCS 
energy). The latter is calculated by using (18), with 
a = (5 = 1, and integrating (21) : 



(E N ~ E sc )bcs 



\ x ±N n (0)Al 



where the extra factor of 1/2 is from the angular inte- 
gral for the d-wave case. Here one would have expected 
a higher gap value than A p (between A p and Ao), but 
the anti-resonant shape of A(E) leads to this simple 
result and is discussed below. The total energy change 
from the normal to SC states is then : 

(En - E sc ) = \N n (0)A 2 p + E v (23) 

If one estimates the new number of gapped pairs as 
N p = N n (0)A p /2y / 2, dividing through by N p gives: 



(E N - E sc ) 1 



N„ 



E v 



This equation matches our expression : Ao = A p • 
with the identification : 



(24) 
-A 



1 A (En — Esc) 



V2 



Neglecting the small imaginary part of Ape, we find 
the energy change (N-^PG) : 



(En — Epq) _ 1 . 



N,, 



V2 



and finally, combining (23) and (24), 
(Epc — E S c) _ E v 



Nn 



(25) 



(26) 



Thus the E v term is precisely the energy gain of the SC 
state with respect to the PG state : the condensation 
energy. The intensive quantity, E v /N p , is of the order 
of k B T c . Since E v is not a 'pairing' energy, we put 
N p = N/2 and obtain : 



A V = 2V2^ 



2V2 k B T c ~ 2.8 k B T c 



in qualitative agreement with the fits. We have thus 
shown that the heuristic parameter A v is consistent 
with a condensation energy, E v , going beyond the BCS 
mean-field value. 

These arguments cannot prove whether A p is a pair- 
ing gap (pre-formed pairs above T c ) or a competing 
order gap (spin density wave, etc.). It is compat- 
ible with many models where states, in the energy 
range ~ A p , are first removed at the Fermi level. 
However, our model suggests that the new quantity, 
A v — 2\/2 E v /N, is not the opening of a second gap 
at the SC transition. Indeed, in a first approximation, 
the total change in KE & PE (i.e. given by the BCS 
Eq. (21)) is identical with the pseudogap contribution : 

(En - Esc)' BC 5 — (^ N — E PG ) (27) 
with the consequence : 

(Epg - EscO'bcs — ( 28 ) 
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However, as previously noted, this is compatible with 
experiment : using STS, the PG width measured at the 
vortex core is about the same magnitude as the SC gap 
between the vortices. The condensation energy term, 
E v — (Epg — Esc) ; is therefore apart, and must be 
due to additional degrees of freedom beyond the mean 
field. 

Behind these results (26-28) is a non-trivial mouve- 
ment of states in the PG — > SC transition [g^. Con- 
sider again Fig. 3 plotting the shape of the gap func- 
tion, and in particular, the role of A„. The states 
at higher energy, E ^ Ed ip , are removed (A(E) ~ 
A p + A v ) while they are added near the anti-resonance 
(A(Eq) ~ A p — A^). Due to a detailed balance of en- 
ergy states in the integrand f(ek) of (21), the effect 
cancels out and we are left with the result Eq. (27). 
This is explicit in Fig. 9 where /(e/c) is plotted as a 
function of €\- (along 8 — 0). We use the previous gap 
expression with a = 1 : 

A(E k ) = A p + (3 A v (I - g(E k )) 

and we compare directly the cases of (3=1, and (3=0. 
As expected, the latter curves cross in such a way as 
the subtended areas are the same. In our view, this 
removal of states at higher energy, also seen in the 
spectral function (Sec. IV), is a significant aspect of 
the problem. A microscopic theory should arrive at a 
clear dependence of E v on this effect. 




20 40 60 BO 100 120 



%(meV) 

Fig. 9. Plot of the integrand, /(efc), of the BCS energy 
equation (21) in units of iV„(0)/(2ir) and for = 0. The 
two curves compare the case of a constant gap, with A p = 
30meV (p = .18), with the case of the gap function, from 
Fig. 3. The latter curve shows a pronounced dip - hump 
feature, but the subtended areas are the same. 

It is possible that, in some instances, the pseudogap 
is larger than the QP peak position: Ape > A p . The 
total energy change, Eq. (23), remains valid but (26) 
becomes : 

(E PG - E SC ) = -jN n (0)(A 2 PG - A 2 p ) + E v 

The change in gap value now opposes the previous 
condensation energy term. (In fact, if Ape is large 
enough, then no transition to the SC state occurs : the 
above vanishes identically.) A disorder potential could 



be such a cause, when amplitude fluctuations are ex- 
pected. Then, the above equation implies a change in 
the number of gapped states at the PG — > SC transi- 
tion and, in this case, E v need not strictly be the order 
parameter. 

The question for further work is to interpret the 
condensation amplitude, A v = 2y/2E v /N. A likely 
candidate is the phase stiffness J relevant to the 
Kosterlitz-Thouless (KT) transition : 

E v /N = (tt/2)J = k B T? T 

For free-electrons in two dimensions, one shows that 
|62| : J = nh 2 /Am, the main point being that J cx n, 
the particle density. From the previous fits to the data, 
we inferred that A v = (p/2p n ) A p , so that it is at 
least proportional to the excess density p. If the KT 
mechanism is involved, then the principal quantity of 
the SC gap function is predicted to be : 

A v = V2ir J = 2V2 k B Tf T 

with a new condition on the phase-stiffness : 



The latter result is satisfactory; the phase rigidity is 
now proportional to both the carrier density and the 
precursor interaction. Put another way, using p' cx 
N n (0) A p as the density of gapped states, we have the 
simple constraint : J cx pp' . Since the latter is valid 
for all p (ranging from to 2po), the dome-shaped SC 
phase diagram is obtained, i.e. pp' cx T c (p). 

In summary, without committing to a specific camp 
in the high-T c problem, the total spectral gap is Ao = 
A p + A^, where A p is the precursor gap and and the 
amplitude A v is proportional to the condensation en- 
ergy. The latter term, and its link to the phase dia- 
gram, points however to the mechanism of 'pre-formed' 
pairs followed by KT condensation. This situation 
of two energy scales was suggested in the work of 
G.Deutscher |79j, where the QP tunneling gap (A p ) 
is distinct from the coherence gap (A v ?) as inferred 
by the Andreev effect. 

Further comments on the gap function 
Consider the question of why the particular QP spec- 
trum, Ek, with the anti-resonant gap function A(-Ek), 
is needed to produce the correct QP-DOS. First, the 
resonance is shown to cause a shift in the quasiparticlc 
states, at the PG to SC transition, from higher energy 
to lower energy, manifested by the sharp QP peaks and 
dip features. Nevertheless, a simple answer for the 
energy change is obtained. All these points indicate 
that the condensate itself is non-trivial [H^] . Further- 
more, the new quantity A v enters the resonance posi- 
tion (Eq ~ Ao), its width and its amplitude (2A V ); a 
property shown to be empirically correct throughout a 
wide doping range. Then it seems unlikely that the res- 
onance is due to the coupling to a quasi-independent 
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collective mode. In the latter case, the width of the 
resonance would reflect the collective mode damping, 
while its amplitude, the strength of the coupling. Here 
these parameters are related to the coherence ampli- 
tude, A v cx Ep/N, possibly tied to the phase-stiffness 
J, and hence to the condensate itself. 

The reason for the energy dependence of the gap 
function is so far unknown from first principles. Aside 
from strong-coupling theory an E dependent gap 
appears in McMillan's proximity effect [80j, two-band 
superconductivity |8 1| . and the asymmetric particle- 
hole case |47| . However, a highly relevant example is 
the detailed work of [82j on the cross-over from BCS 
to Bosc-Einstein condensation; but the gap has no res- 
onance effect there. 

The detailed form of the gap function, or self-energy 
as in Section IV, represents a novel pairing interac- 
tion. It is possible that there is a strong feedback effect 
between the precursor state, with gap Ape, and the 
freezing of phase fluctuations, of amplitude A v . Then 
there is a complete renormalization of the gap function 
A(Ek), where Ek is the 'exact' quasiparticle spectrum 
(final state). The effect is strong enough so that if 
Ek ~ Ao, the pairing energy is reduced, A ~ A p — A v , 
compared to the low-energy value A ~ A p . On the 
contrary, at higher energy, the interaction is strong 
A ~ A p + A^. This type of variation is at the heart 
of the QP-DOS shape. 



which has the expected Lorentzian form with a peak 
at ReEk = E and half- width ImE^. Such an expres- 
sion neglects many- body effects that are important at 
higher energy 0,^1- 
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IV. SPECTRAL FUNCTION 

A different, yet equivalent, view of the same prob- 
lem is given by the spectral function, A(k,E). It is 
more fundamental since there is no sum over momenta, 
hence an ideal quasiparticle with wave-vector k is di- 
rectly a Dirac peak at the position E = E^. It is also 
the fundamental quantity describing the ARPES mea- 
surement. 

Here we discuss A(k,E) as given by equation (8), 
with the quasiparticle dispersion (7), the total gap 
function : 

A(E k ) = aA PG (1 - 0) + A S c(£k) 

and with Asc(Ek) the usual resonant form given in 
(10). The gap broadening parameters are as in Sec. 
II, i.e. A p ->■ A p (1 - iS) and A PG ->■ A p (1 - iS') in 
Asc and Ape, respectively. The final broadening pa- 
rameter is r Dynes, and all three have been determined 
by the fits in Sec. II. 

Since A(k, E) is calculated along the real axis, it 
follows that our dispersion law E k = -Ek(efe), Eq. (7), 
is complex. Considering the anti-nodal direction, the 
simple model gives : 



A(k,E) = 



Im Ek 



1 



7T (E-ReE k y 



ImEl 



(29) 



Fig. 10. Direct comparison of the spectral functions of 
the PG state (upper panel) and of the SC state (lower 
panel), viewed as density plots in the (E,ek) plane. The 
gap parameters correspond to p = .18 : A p = 30meV and 
A v — 17.8 meV. Other relevant parameters are: S = .06, 
5' = .2 and T = .06 A p . 

The dispersion law, ReEk, which is hyperbolic in the PG 
case, has a distinctive kink (at D) in the SC case, causing 
the dip feature in the DOS. States are removed at high 
energy from the PG spectral density (compare M and M'), 
and are added at low energy (between P and D), in the SC 
spectrum. 

The spectral function (29) is presented in grey scale 
in Fig. 10 in the (E, e&) plane and for the case where 
p = .18. Note that the gap function is an explicit 
function of the QP energy, E^, and not ek', thus one 
must first find numerically the (complex) dispersion 
law Ek(ek) in order to apply (29). On the lower panel 
of Fig. 10 we give the spectral function for the SC state 
and compare directly with the one for the PG state 
(upper panel). For the SC state, the visibly sharp 
peak position describes a continuous curve that de- 
viates significantly from the BCS hyperbolic law. The 
hyperbolic case is seen in the upper panel, since it is 
inherent in the model for the pseudogap. For small e 
both curves begin near the gap value E ~ A p = 30 
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meV, as expected. The QP peak position (SC state) 
begins (at P) with a slightly smaller slope than the 
BCS case and nearly reaches the E = e k diagonal. 
This is followed by a sharp kink (near efc = 60 meV, 
at D) before joining the expected asymptote for a gap 
value of Ao ~ 48meV. All of these effects are totally 
absent in the BCS framework. 

The kink feature in the dispersion law corresponds 
to the dip in the QP-DOS. One can read its position, 
D, either on the E axis or on the efc axis: in the for- 
mer, Edi P ~ A p + 2 A v ~ 64 meV, while in the latter 
£fc ~ 2y / AoA^ ~ 60 meV (taking the point of steep- 
est slope). The anti-resonance in the gap function is 
revealed here by the near touching of the diagonal by 
the dispersion law at the energy E ~ Ao ~ 48meV. 
(At this point efc ~ E.) Concerning the peak intensity, 
the sharpest QP peaks are not at efc = 0, E = A p (P in 
Fig. 10) but at slightly higher energy, near 40meV (at 
P'). At the kink, the peak intensity drops significantly 
(highest broadening), giving a quite different picture 
of the peak to dip feature. 



O 
Q 



i 




Energy (meV) 

Fig. 11. Tunneling DOS (along 8 = 0) calculated using 
the exact formula (9), solid line, and the discrete sum over 
efc of the spectral density of Fig. 10, triangular dots. 

Although ARPES can only access the occupied 
states, contrary to tunneling spectroscopy, the two 
techniques probe the same spectral function (Fig. 10). 
Indeed, integrating A(k, E) over efc, for a given E, 
yields the partial tunneling DOS, n s (E,0), as given 
by Eq. (6). In ARPES one can fix k, hence e^, and 
measure A(k, E) as a function of E (energy distribu- 
tion curve). Moreover, the dispersion of the QP states 
with k|| can be probed. Although these two techniques 
have different spatial and energy resolutions, it would 
be of high interest to establish the common spectral 
function that matches both the ARPES and the tun- 
neling data. 

As a simple check of the tunneling DOS, in Fig. 11 
we plot both the exact partial DOS (for 9 = 0), using 
eq.(9), and the direct sum over states of the spectral 
density of Fig. 9. Given the crude integration method 
in the latter case, the agreement is satisfactory. All 
experimentally observed features are found : the pro- 
nounced QP peak and the dip, as well as the states 
within the gap from the Dynes broadening. 

The upper panel of Fig. 10 shows the spectral func- 



tion corresponding to the PG state, or (3 = in the 
gap function, for direct comparison. Notice the sig- 
nificant smearing near e& ~ 0, or E ~ A p : as ex- 
pected, there are no well-defined quasiparticle peaks. 
The overall dispersion shape is still hyperbolic and fol- 
lows Ek " 



e 2 + A 2 as imposed by the model. The 
asymptotic line is closer to the diagonal as compared 
to the SC case, where the total high-energy spectral 
gap is larger, i.e. Ek ~ \J e 2 + A 2 , (compare M' with 
M in Fig. 10). This gives a direct view of the mouve- 
ment of the states in the PG — > SC transition : states 
well above the PG gap, for Ek > Edi P (upper panel), 
are removed and a high density of quasiparticle states 
exist, for Ek < Edi P , in the SC state (lower panel). 

Quasiparticle self- energy 

The final objective of this work is to show that the E- 
dependant gap function is equivalent to a self-energy. 
In the SC state, the carriers of energy efc are modified 
by the presence of Ek(efc) : 

A(k, E) = \ Im ^ , ~ (30) 

tt E - efc - 2jk(£fe) 

Thus, setting S = gives the normal state spectrum, 
while for ideal BCS quasiparticles it is with Ek(efc) = 
^/e 2 . + A 2 — efc. Finally, the Dynes lifetime broadening 
is obtained if £ has an imaginary part (Im £ = T). 

For the general case, where the SC state is given 
by the gap function, comparing (30) with our original 
expression (8), gives: 



^k(efc) — E k (e k ) — e fc 



(31) 



However, the final state QP energy, i?k(efe), is a solu- 
tion of the implicit equation E^ = v/e? + Ak(-Ek) 2 + 
iT, and thus an analytical expression for £k(efc) can- 
not be given. The computation of -Ek(efc) is quite 
straightforward, leading to the complex valued £k(efc) 
and hence to the spectral function (30). Formally, the 
quasiparticle resonance position is given by ReE^ = 
efc + i?e£k(e/c), and its width by ImE^ = 7r7iSk(efe)- 
Note that using the definition of E^, and the simple 
expression (31), yields: 



Sk(efc) = iT 



A k (£ k 



Ek - e k + iT 



(32) 



which is similar to the form considered by 0, 
Thus, the independent gap function and the self- 
energy function are equivalent. 

The self-energy is plotted in Fig. 12a using the same 
parameters as those of Figs. 10 and 11. The darker 
line corresponds to the SC case using the resonant gap 
function, the gray line is for the PG model, with a 
constant gap (~ A p ). For small efc both curves begin 
at the value: i?e£(0) ~ A p ~ 30 meV, as expected. 
However, the SC curve descends sharply towards zero, 
as compared to the PG one, and reaches a minimum 
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near e k ~ 2^J A p A^,. This corresponds to the anti- 
resonance position in the gap function (E = Aq). The 
self-energy then rises abruptly, crosses the PG refer- 
ence curve and, with a further change in concavity, 
decays as a function The dip position in the QP- 
DOS is found by the maximum slope of Re £(£&); it is 
near the crossing point at e k — 2^/Ao A v . Note that 
this dip position coincides with the highest broadening, 
as revealed by ImS(efc) in the lower panel of Fig. 12a. 
The lowest broadening is at t k ~ 40 meV, giving the 
highest QP peaks in the spectral function, and the wide 
peaks seen in the DOS. 



(a) (b) 




20 40 60 SO 100 120 -150 -100 -50 50 100 150 



Fig. 12a. Complex self-energy plotted against e k along 
8 — 0. Dark line - SC self-energy; gray line - PG case. 

The strong dip, where ReE(tfc) nearly vanishes, is fol- 
lowed by a crossing-point with the PG line, then slowly 
decays at higher energy. The sudden rise in ReE(efc), at 
a maximum of ImE(ejb), is at the energy e k — 2y/ Ao A v , 
and causes the DOS dip. The sharpest QP peaks are not 
at e k ~ 0, E ~ A p , but nearer ~ 40 meV. 

Fig. 12b. Plot of the pair amplitude 2 F k = A k (E)/E k , 
using the gap function (18), and compared to the BCS 
value, along the anti-nodal direction. (All parameters are 
the same as in Figs. 8-10.) 

In conclusion, the self-energy of Fig. 12a is yet an- 
other way of comparing the change of states from the 
PG to SC cases. Indeed, the self-energy of the SC 
state is smaller than the PG one for small energies, 
up to the crossing point, when the relative magni- 
tude inverts. Our previous observation that states 
are removed at large energy from the PG state, for 
£/t 3> 2-y/Ao A v , while they are added to the SC state, 
for £fe < 2-^/Ao Ay, is seen here in a different way. The 
strong minimum of i?eS, near e k ~ 2y/ A p A v , also re- 
veals the condensation amplitude. The sharpness of 
these features is an indication that ARPES may be a 
better technique to extract these two parameters, as 
compared to quasiparticle tunneling. 

In Fig. 12b we plot the pair amplitude Fk = A k /2E k 
as a function of |Z£j, taking the real part. The case 
where the SC gap function is used is directly compared 
to the conventional pair amplitude, with a constant 
gap. The obvious dip features, near ± 50 meV, are 
due to the same minimum discussed previously. The 
investigation of this function is yet to be done in the 
case of cuprates; it requires the Andreev effect or pair 



tunneling, and not quasiparticle tunneling, in order to 
be probed. 

V. CONCLUSIONS 

Three principal questions have been addressed : does 
the tunneling DOS reveal the order parameter and if 
so, in what manner, and finally, does it connect to the 
phase diagram? We have tried to answer them from 
the heuristic point of view, i.e. by using our resonant 
gap function and examining the available data. This 
-E-dependant gap produces the well-known quasipar- 
ticle peak shapes, followed by the dip structures at 
higher energies. No Fermi surface anisotropy, no van 
Hove singularity, no strong-coupling are needed; only 
simple d-wave quasiparticles, with a modified spec- 
trum, are considered. 

The model contains two fundamental energy scales, 
whose values are determined by the fits to the tunnel- 
ing spectra: A p , the quasiparticle peak position, and 
a second energy, A v . We found the resonance (along 
the anti-nodal direction) to be at the fixed position 
Ao = A p + A v , coinciding with the total spectral gap 
for high energies, i.e. E k — \/ e\. + A'q. Thus, the de- 
scription of the SC state depends uniquely on this new 
energy scale : A v . 

To go a step further, we showed that these two en- 
ergy scales fit the phase diagram as a function of car- 
rier density, p. While A p decreases linearly with dop- 
ing, A v follows T c , with A v — 2.3fcsT c , and merges 
smoothly with A p at high doping. Thus the gap func- 
tion scales perfectly with the carrier density; one can 
then generate the QP-DOS shape for arbitrary p. The 
two energy scales have a simple link : A v cx p A p , valid 
in the range of the available data. If there is a criti- 
cal point, it would be at maximal doping where both 
energy scales vanish. 

The role of A v as a new condensation amplitude is 
compelling. To check this, we estimated the energy 
changes associated with the three states (normal, PG, 
SC). A simple interpretation for each term of the con- 
straint A = A p + A v is, respectively : the total energy 
gain (per pair), the energy gain due to the precursor 
gap (PG state) and the condensation energy per pair 
E v /N p . However, the number of gapped pairs turns 
out to be N p oc A p /2 and not N p oc Aq/2, which is 
thus in support of the preformed-pair model. For a 
competing order gap, N p is to be simply taken as half 
the number of gapped states. In either case our new 
amplitude turns out to be A v = ^/2E V /N P . It is there- 
fore justified, for both scenarios, that this parameter 
is proportional to T c , as found independently from the 
fits to the data. 

Since the mouvement of states in the PG to SC tran- 
sition remains subtle, i.e. not the direct opening of a 
second (SC) gap, we took a closer look at the disper- 
sion lawi?k(efc)- We discussed the spectral function, 
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and the corresponding self-energy, where the change 
of the quasiparticle states can be seen more readily 
than in the DOS. Indeed, it allows to identify features 
of the QP spectrum as a function of both E and e^. 
One finds that in the PG state quasiparticles of higher 
energy, tk S> 2WAo A^, end up as quasiparticle peaks 
at lower energy, < 2-^/Ao A v , in the SC state. The 
cross-over is precisely the characteristic energy where 
the dip feature is seen, either in the self-energy or in 
the gap function. By a different analysis than of the 
DOS, we found that this mouvement of the states de- 
pends uniquely on A v . 

Our model shows that the principal change of kinetic 
and potential energies (as in BCS) is associated with 
the opening of the precursor gap Ape — A p , at the 
Fermi level (normal state to PG state transition). At 
the SC transition, the quantity E v is therefore linked 
to other degrees of freedom not accounted for in the 
BCS model. This possibility is accentuated by the res- 
onant character of the gap function. While a particular 
collective effect is perhaps responsible, a likely candi- 
date for the condensation energy is the phase stiffness, 
J, in the Kosterlitz-Thouless transition. In this con- 
text, we get an immediate relation between our param- 
eter and J : A v = ^/2n J which is then proportional to 
the KT transition temperature. A relevant prediction 
is thus: J cx pA p , should the KT mechanism fit the 
phase diagram, for a wide range of doping. 

The three questions above have therefore been given 
a tentative answer. One could raise some objection 
to introducing an ad hoc gap function, where a micro- 
scopic theory is missing. Some speculations concerning 
the physical origin of the gap function have been given 
in the closing paragraphs of Sec. III. Our main goal 
of matching the observed experimental QP-DOS, and 
their characteristics as a function of the carrier den- 
sity, works without resorting to additional extraneous 
factors. Moreover, the fits can now be done with a 
minimal set of parameters. Finally, the two key pa- 
rameters, the QP peak position A p , and the coherence 
amplitude A^,, have been connected to observed phys- 
ical quantities. 

One of the authors (W.S.) thanks Amir Kohen & 
Th. Proslier, for interesting and helpful discussions. 
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